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This  paper  introduced  a  one-dimensional  analytical  model  to  quantitatively  examine  the  humidifying 
capacity  of  a  Nation™  membrane  humidifier.  Water  permeability  for  the  Nation™  membrane  was  calcu¬ 
lated,  depending  on  the  thickness  of  the  Nation™  membrane.  Relative  humidity  for  a  carrier  gas  passed 
through  the  Nation™  membrane  humidifier  was  simulated  at  various  gas  flow  rates,  which  showed  good 
agreement  with  experimental  data.  Simulations  were  also  conducted  to  predict  humidity  levels  as  a  func¬ 
tion  of  length  and  height  of  the  gas  flow  channel.  This  analytical  model  can  be  used  in  future  work  to 
properly  design  Nation™  membrane  humidifiers  for  use  in  PEM  fuel  cell  systems. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proton  exchange  membrane  (PEM)  fuel  cells  are  highly  efficient 
energy-conversion  devices  for  automotive  and  stationary  applica¬ 
tions.  Typically,  the  humidified  fuel  and  oxidant  are  introduced  to 
the  anode  and  the  cathode  compartments  in  order  to  make  the 
membrane  fully  hydrated  in  a  wide  current  range  [1-3]. 

Past  research  has  focused  on  developing  various  humidifying 
methods  to  effectively  supply  humid  reactants  into  the  PEM  fuel 
cells  [4-9].  Reid  [4]  attempted  to  humidify  a  reactant  by  passing 
it  through  a  heated  water  reservoir,  while  Nguyen  and  White  [5] 
suggested  direct  water  spraying  into  the  gas  stream.  Chow  and 
Wozniczka  [6]  developed  a  membrane  humidifier  capable  of  serv¬ 
ing  water  vapor  to  a  reactant.  As  illustrated  in  Fig.  1 ,  the  membrane 
is  sandwiched  between  the  gas  flow  channel  (i.e.,  permeate  side) 
and  the  water  (or  humid  gas)  flow  channel  (i.e.,  feed  side).  Water 
permeates  through  the  membrane  and  evaporates  into  the  gas 
stream  due  to  a  lower  partial  pressure  exerted  on  the  permeate 
side  of  the  membrane  [10]. 

Recent  simulation  studies  [11-14]  yielded  a  deeper  understand¬ 
ing  of  how  water  moves  into  a  carrier  gas  through  the  Nation™ 
membrane  for  various  operating  conditions:  gas  flow  rate,  gas  or 
water  temperature,  water  pressure.  The  models  in  these  studies 
were  proposed  to  calculate  water  content  in  a  carrier  gas  at  the 
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outlet,  in  terms  of  Nation™  membrane  humidifier  performance. 
Although  these  models  contribute  to  the  analysis  and  prediction 
of  water  transport  through  the  Nation™  membrane  in  a  humidi¬ 
fier,  they  are  difficult  to  apply  in  practice  because  of  their  complex 
forms.  Hence,  a  simplified  model  is  necessary  and  helpful  to  evalu¬ 
ate  humidification  using  the  Nation™  membrane  for  PEM  fuel  cell 
applications.  In  this  study,  a  one-dimensional  analytical  expression 
was  developed.  Water  permeability  was  calculated  based  on  pre¬ 
vious  experimental  results  [9].  The  relative  humidity  for  a  carrier 
gas  passed  through  the  Nation™  membrane  humidifier  at  various 
flow  rates  was  predicted  using  an  analytical  model  and  validated  by 
comparing  with  experimental  data.  The  effects  of  the  dimensions 
of  the  gas  flow  channel  on  the  humidifying  capacity  of  the  Nation™ 
membrane  humidifier  were  also  examined. 

2.  Mathematical  model 

The  mass  flux  of  water  from  the  water  flow  channel  to  the  gas 
flow  channel  through  the  membrane  in  Fig.  1  is  expressed  as  [10,15] 

N™em  —  kw(*w?WPw,0  -yWP)  (i) 

where  /<w  is  the  water  permeability,  xw  is  the  mole  fraction  of  the 
water  in  liquid  phase,  yw  is  the  activity  coefficient  of  the  water, 
pW) o  is  the  saturated  pressure  of  the  water  at  a  given  temperature, 
yw  is  the  mole  fraction  of  the  water  in  gas  phase,  and  P  is  the  total 
pressure.  Based  on  the  assumptions  that  (i)  the  gases  behave  as 
ideal  ones,  (ii)  the  membrane  humidifier  is  isothermal,  (iii)  the 
total  pressure  remains  constant  (i.e.,  ambient  pressure),  (iv)  the 
partial  pressure  of  water  vapor  changes  only  in  the  z-direction,  (v) 
the  water  permeability  is  constant,  and  (vi)  the  gas  flow  channel 
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is  straight,  the  mass  balance  on  water  vapor  along  the  gas  flow 
channel  can  be  established  as  follows: 


Nw\WH  -  Nw|  z+AzWH  +  WPw,o-Pw)WAz  =  0  (2) 


where  W  is  the  width  of  the  gas  flow  channel,  H  is  the  height  of  the 
gas  flow  channel,  and  A z  is  the  thickness  of  an  intersection  of  the 
gas  flow  channel.  Dividing  Eq.  (2)  by  WHAz  and  taking  the  limit  as 
A z  approaches  zero  gives 


dNw 

dz 


/<w 

IT 


(Pw,0  ~  P w) 


0) 


The  partial  pressure  of  water  vapor  is  a  function  of  both  x  and  z. 
However,  the  partial  pressure  of  water  vapor  can  be  considered  vir¬ 
tually  as  a  function  of  z  alone  since  the  velocity  profile  of  a  carrier  gas 
is  not  significantly  affected  by  the  evaporation  across  cross-section 
of  the  gas  flow  channel  [16].  When  the  flow  in  the  gas  channel 
is  dominated  by  convection  (i.e.,  Pe  =  vL/Dw-g»  1),  the  constitute 
equation  for  water  flux  in  the  gas  flow  channel  is 


A/w  —  — CDW-g 


dXy\j 

dz 


+  Xw(iVw  +  Ng) 


CwV 


(4) 


with  cw  = 


(5) 


where  Dw-g  is  the  water/gas  binary  diffusion  coefficient,  Ng  is  the 
gas  flux,  cw  is  the  water  concentration,  v  is  the  gas  flow  rate,  R  is  the 
gas  constant,  and  T is  the  temperature.  Substituting  Eqs.  (4)  and  (5) 
into  Eq.  (3)  yields 


dpw 

dz 


kwRT 

vH 


(Pw,0  Pw) 


(6) 


If  the  partial  pressure  of  the  water  at  the  gas  flow  channel  inlet 
Pwlet  is  known,  Eq.  (6)  can  be  solved  analytically  as  follows: 


Pw 

Pw,0 


=  1  exp 


3.  Results  and  discussion 


3  A.  Determination  of  water  permeability 

The  relative  humidity  Hr  is  defined  as 

Hr  =  —  X  100  (8) 

Pw,0 

As  shown  in  Fig.  2,  the  relative  humidity  Hr  measured  at  various 
air  flow  rates  was  plotted  in  the  plane  of  -ln(l  -pw/pw,o)  versus 
v-1,  to  calculate  water  permeability  kw  for  a  Nation™  membrane 
with  different  thicknesses.  The  experiment  performed  to  obtain  the 


Fig.  2.  Determination  of  the  water  permeability  for  Nation™  membranes,  (inset) 
The  water  permeability  coefficient  for  a  Nafion™  membrane. 


relative  humidity  data  in  Fig.  2  is  fully  described  in  [9].  In  Fig.  2,  the 
value  of  kw  for  the  Nafion™  membrane  is  determined  by  modifying 
Eq.  (7).  i.e., 


■In  f  1  -  ZiA  = 

V  Pw,0  J  vH 


O) 


where  L  is  the  length  of  gas  flow  channel  in  the  Nafion™  mem- 
brane  humidifier  (I  =  0.297  m).  The  values  of  kw  calculated  from  Eq. 
(9)  are  8.9  x  10-5,  5.3  x  10-5,  and  2.3  x  10-5  kgm-2  s-1  atm-1  for 
Nafion™  112, 115,  and  117,  respectively.  These  results  indicate  that 
values  of  kw  for  the  Nafion™  membrane  are  strongly  related  to  the 
thickness. 

The  process  of  membrane  humidification  includes  the  diffu¬ 
sion  of  water  through  the  membrane  and  interfacial  mass-transfer, 
which  includes  sorption  into  the  membrane  on  the  feed  side  and 
desorption  into  a  carrier  gas  on  the  permeate  side  [10,11,17].  If  the 
interfacial  mass-transfer  is  dominant  in  the  membrane  humidifica¬ 
tion,  the  water  permeability  is  highly  related  to  gas  flow  rate.  This  is 
mainly  because  the  thickness  of  the  boundary  layer  at  the  permeate 
side  (i.e.,  the  interface  between  a  carrier  gas  and  the  membrane)  is 
largely  affected  by  gas  flow  rate,  when  the  boundary  layer  at  the 
feed  side  is  assumed  to  be  negligibly  thin.  On  the  other  hand,  if  the 
diffusion  through  the  membrane  is  dominant  in  humidification,  the 
value  of  kw  is  inversely  proportional  to  thickness  of  the  membrane. 

The  permeability  coefficient  of  water  across  a  membrane  Pw  is 
the  product  of  water  permeability  kw  and  thickness  of  a  mem¬ 
brane  8mem  [10,15].  As  seen  in  Fig.  2,  the  values  of  Pw  for  a 
Nafion™  membrane  with  different  thicknesses  are  between  4.0 
and  6.6  x  10-9  kgm-1  s-1  atm-1  (see  the  inset).  The  inconsistency 
of  Pw  values  may  originate  primarily  from  the  theoretical  assump¬ 
tion  that  the  mass  flux  of  water  vapor  is  expressed  as  kw(pw,  o  -  Pw) 
where  kw  is  not  a  function  of  position  but  constant  and  pw  is  a  func¬ 
tion  ofz  alone.  This  issue  will  be  addressed  in  future  work.  However, 
even  if  the  values  of  Pw  vary  somewhat,  it  is  clear  that  diffusion  is 
the  dominant  process  in  determining  kw,  since  kw  is  not  appreciably 
affected  by  gas  flow  rate. 


3.2.  Effect  of  gas  flow  rate 

Fig.  3  displays  the  values  of  Hr  for  a  carrier  gas  passed  through 
a  Nafion™  membrane  humidifier  at  various  flow  rates  using  Eq. 
(7).  The  simulations  were  performed  at  60  °C.  The  values  of  kw 
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Fig.  3.  Relative  humidity  for  a  carrier  gas  passed  through  a  Nation™  membrane 
humidifier  as  a  function  of  gas  flow  rate. 

for  Nafion™  111  (=25  fxm)  and  1135  (=87  |jim)  were  calculated  as 
2.0  x  10-4  and  6.7  x  10-5  kgm-2  s-1  atm-1,  which  were  calculated 
using  Pw  =  5.0  x  10-9  kgm-1  s-1  atm-1,  and  the  values  of  L  and  H 
were  0.297  m  and  1  mm.  Our  previous  experimental  data  are  also 
shown  in  Fig.  3  [9].  The  simulation  results  shown  in  the  figure  agree 
with  the  values  of  HR  measured  at  various  air  flow  rates.  When  the 
gas  flow  rate  v  increases,  HR  values  for  a  carrier  gas  decay  drasti¬ 
cally  at  low  flow  rates,  especially  for  a  thicker  Nafion™  membrane. 
Hr  values  at  v  =  0.333  m  s-1  were  99.2%  for  Nafion™  111,  88.6%  for 
Nafion™  112,  80.2%  for  Nafion™  1135, 72.3%  for  Nafion™  115,  and 
43.0%  for  Nafion™  117,  respectively.  At  high  flow  rates,  HR  decreases 
slowly,  showing  smaller  differences  between  Nafion™  membranes. 
These  results  imply  that  a  carrier  gas  is  not  sufficiently  humidified 
at  high  flow  rates  due  to  short  residence  time  in  the  membrane 
humidifier,  even  if  a  thinner  Nafion™  membrane  is  used. 

3.3.  Effect  of  gas  flow  channel  length 

In  a  membrane  humidifier,  a  carrier  gas  supplied  into  the  gas 
flow  channel  is  humidified  when  it  moves  in  the  z-direction,  as 
portrayed  in  Fig.  1 .  Hence,  it  is  very  important  to  predict  the  humid¬ 
ifying  capacity  of  Nafion™  membranes  as  a  function  of  length  of  the 
gas  flow  channel,  which  is  coupled  with  the  volume  of  the  Nafion™ 
membrane  humidifier.  Fig.  4(a)  delineates  the  values  of  HR  along  the 
gas  flow  channel  for  a  carrier  gas  passed  through  Nafion™  mem¬ 
brane  humidifier.  The  simulations  were  conducted  at  60  °C,  and  the 
values  of  v  and  H  were  0.333  m  s-1  and  1  mm.  As  represented  in 
Fig.  4(a),  the  lengths  of  the  gas  flow  channel  necessary  to  obtain  HR 
of  80%  for  a  carrier  gas  were  0.10,  0.22,  0.29,  0.37,  and  0.85  m  for 
Nafion™  111,  112, 1135, 115,  and  117,  respectively. 

A  carrier  gas  introduced  into  a  Nafion™  membrane  humidifier 
may  contain  water  vapor,  relying  on  methods  such  as  the  electrol¬ 
ysis  and  steam  reforming  process  to  produce  reactants.  Fig.  4(b) 
specifies  the  values  of  HR  along  the  gas  flow  channel  using  differ¬ 
ent  Nafion™  membranes  at  HR  inlet  of  40%.  Comparing  the  profiles 
of  Hr  for  different  Nafion™  membranes  in  Fig.  4(b),  the  pattern 
of  Hr  values  is  similar  to  that  shown  in  Fig.  4(a).  However,  the  gas 
flow  channel  lengths  to  reach  HR  of  80%  are  shorter  by  ca.  60%  than 
those  with  vapor-free  gas  at  the  inlet  (i.e.,  Fig.  4(a)).  Fig.  5  depicts 
the  values  of  HR  as  a  function  of  length  of  the  gas  flow  channel  in 
Nafion™  112  membrane  humidifier,  when  the  values  of  HR  at  the 
inlet  are  0, 20, 40,  60,  and  80%,  respectively.  Note  that  the  gradients 


Fig.  4.  Relative  humidity  for  a  carrier  gas  passed  through  a  Nafion™  membrane 
humidifier  as  a  function  of  gas  flow  channel  length  at  (a)  HRjinIet  =  0%  and  (b) 
HR,inlet  =  40%. 


Fig.  5.  Relative  humidity  for  a  carrier  gas  along  the  gas  flow  channel  at  various  HR  inlet 
for  Nafion™  112. 
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Table  1 

Maximum  gas  flow  rates  to  obtain  Hr  of  80%  using  Nafion™  membrane  humidifiers 
and  the  currents  generated  under  stoichiometry  with  A,air  =  3  (operating  temperature 
in  Nafion™  membrane  humidifiers  and  a  PEM  fuel  cell  is  60  °C). 


Fig.  6.  Relative  humidity  for  a  carrier  gas  passed  through  a  Nafion™  membrane 
humidifier  as  a  function  of  gas  flow  channel  height. 


heights  using  Eq.  (7).  The  simulations  were  performed  at  60  °C,  and 
the  values  of  v  and  L  were  0.333 ms-1  and  0.297m,  respectively. 
The  value  of  Hr  at  the  outlet  of  the  Nafion™  membrane  humid¬ 
ifier  remains  near  100%  and  then  decreases  significantly  as  gas 
flow  channel  height  increases,  depending  on  the  thickness  of  the 
Nafion™  membrane.  Interestingly,  the  value  of  HR  for  a  carrier  gas 
passed  through  the  Nafion™  111  membrane  humidifier  is  higher 
than  99%,  when  the  gas  flow  channel  height  ranges  from  0  to  1  mm. 
As  illustrated  in  Fig.  6,  when  the  thickness  of  Nafion™  membrane 
varies  between  25  and  175  p,m,  the  value  of  Hr  from  the  Nafion™ 
membrane  humidifier  with  H  =  2  mm  drops  from  91.3%  to  24.5%.  The 
simulations  in  Fig.  4(a)  and  6  clarify  that  in  the  case  of  Nafion™  111, 
the  length  and  height  of  the  gas  flow  channel  necessary  to  reach  Hr 
of  80%  for  a  carrier  gas  at  v  =  0.333  m  s-1  are  ca.  8.5  times  smaller 
and  larger  than  those  in  the  case  of  Nafion™  117,  respectively.  That 
is,  the  minimum  volume  of  the  gas  flow  channel  for  Nafion™  111 
to  operate  a  fuel  cell  at  a  given  current  without  a  significant  drop 
in  ionic  conductivity  of  the  membrane  is  ca.  8.5  times  smaller  than 
that  for  Nafion™  117. 

4.  Conclusion 

The  results  of  this  study  indicated  that  a  thinner  Nafion™  mem¬ 
brane  leads  to  higher  relative  humidity  for  a  carrier  gas,  since  water 
permeability  depends  heavily  on  the  thickness  of  a  Nafion™  mem¬ 
brane.  Relative  humidity  for  a  carrier  gas  passed  through  a  Nafion™ 
membrane  humidifier  was  calculated  using  an  analytical  expres¬ 
sion  along  the  gas  flow  channel,  as  a  function  of  gas  flow  rate,  length, 
and  height  of  the  gas  flow  channel.  This  analytical  model  can  be  uti¬ 
lized  to  effectively  design  Nafion™  membrane  humidifier  for  PEM 
fuel  cells. 
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of  Hr  near  the  inlet  decreases  with  increasing  HR  inlet,  whereas  the 
lengths  of  the  gas  flow  channel  needed  to  reach  Hr  of  100%  decrease 
with  higher  HR  inIet. 

According  to  Wilkinson  et  al.  [18],  relative  humidity  of  80%  for 
an  oxidant  introduced  to  the  cathode  inlet  is  high  enough  to  suf¬ 
ficiently  maintain  ionic  conductivity  of  the  membrane.  Assuming 
that  the  stoichiometry  of  the  fully  humidified  hydrogen  is  suffi¬ 
ciently  high  at  the  anode  compartment  in  a  PEM  fuel  cell  (Ah2  »  1 ), 
the  maximum  gas  flow  rate  vmax  to  reach  Hr  of  80%  for  a  carrier  gas 
using  Nafion™  membrane  humidifier  (L  =  0.297  m)  and  the  current 
generated  from  a  PEM  fuel  cell  under  constant  stoichiometry  with 
Aair  =  3  are  summarized  in  Table  1.  When  air  is  used  as  an  oxidant, 
the  current  generated  under  constant  stoichiometry  is  calculated 
as  follows  [19] 


_  4pmaxAchEP 

i^air 


(10) 


where  Ach  is  the  cross-section  area  of  the  gas  flow  channel  at  the 
cathode  compartment,  F  is  the  Faraday  constant,  P  is  the  total  pres¬ 
sure  in  a  PEM  fuel  cell,  respectively.  As  listed  in  Table  1,  the  value  of 
/air  for  Nafion™  111  is  8.7  times  higher  than  that  for  Nafion™  117. 
Using  Eqs.  (7)  and  (10),  the  minimum  length  of  the  gas  flow  channel 
in  the  Nafion™  membrane  humidifier  for  demanding  maximum 
current  can  be  quantitatively  estimated. 


3.4.  Effect  of  gas  flow  channel  height 
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Fig.  6  exhibits  the  values  of  Hr  for  a  carrier  gas  passed  through 
a  Nafion™  membrane  humidifier  with  different  gas  flow  channel 


